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Abstract 

The s-point correlation function of a Gaussian Hermitian random matrix 
theory, with an external source tuned to generate a multi-critical singular¬ 
ity, provides the intersection numbers of the moduli space for the p-th spin 
curves through a duality identity. For one marked point, the intersection 
numbers are expressed to all order in the genus by Bessel functions. The 
matrix models for the Lie algebras of 0{N) and Sp{N) provide the intersec¬ 
tion numbers of non-orientable surfaces. The Kontsevich-Penner model, and 
higher p-th Airy matrix model with a logarithmic potential, are investigated 
for the open intersection numbers, which describe the topological invariants 
of non-orientable surfaces with boundaries. String equations for open/closed 
Riemann surface are derived from the structure of the s-point correlation 
functions. The Gromov-Witten invariants of CP^ model are evaluated for 
one marked point as an application of the present method. 
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1 Introduction 


The intersection numbers of moduli space of Riemann surface are topologi¬ 
cal invariants, which are closely related to universal singularities in random 
matrix models with an external source. By tuning appropriately the matrix 
source, multi-critical behaviors are obtained at the edge of the density of 
state [1, 2, 3]. They are described by Airy and higher Airy kernels for the 
p-th degenerate singularity. The correlation functions for the p-th singular¬ 
ity turn out to be generating functions for the intersection numbers of p-spin 
curves on Riemann surfaces [4, 5, 6, 7]. In this article, we extend our previous 
work on the calculation of the intersection numbers of p spin curves. This 
technique relies on a duality [4, 6] which is specihed in the next section. 

Recently, the open intersection theory with boundaries has been investi¬ 
gated in [16, 17, 18]. The generating matrix model for these open intersec¬ 
tion numbers is the Kontsevich-Penner model [21, 22], namely Kontsevich’s 
model with a logarithmic potential, studied before in [19, 20]. The Virasoro 
equations for this case provide a different structure from KdV, since the Rie¬ 
mann surface with boundaries is no longer orientable. The expansion of the 
logarithmic matrix model, like a ribbon graph expansion of Kontsevich ma¬ 
trix model, provides the non-orientable surface. The situation is similar to 
that of non-orientable surfaces generated by the matrix models of 0{N) and 
Sp{N) Lie algebras [12]. We compute here the open intersection numbers, 
and higher p-th spin curves, including such Lie algebras. 

The composition of this article is the following. In section 2, several du¬ 
alities exchanging the size of the random matrices in an external matrix with 
the number of points in correlation functions, are recalled. The dehnition 
of the intersection numbers is briefly recalled. In section 3, the intersec¬ 
tion numbers for p = 2, 3,4, 5, and p = —1 (Euler characteristics) and one 
marked point, are explicitly given in terms of Bessel functions. In section 4, 
the intersection numbers for the Lie algebras of 0(2N), 0(2N-|-1), and Sp(N) 
are discussed. The Euler characteristics for those non-orientable surfaces 
follow. In section 5, the open intersection numbers are computed for one 
marked point from the Kontsevich-Penner model , and the relation to the 
0{N) model for non-orientable intersection numbers is discussed. In section 
6, the intersection numbers for multiple marked points are evaluated and 
the relation to Virasoro equations is examined. Section 7 is application of 
the present method to the Gromov-Witten invariants of CP^. Section 8 is 
devoted to discussions. In an appendix we study the case of p-spin curves 
for open Riemann surfaces. 
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2 Dualities 


• GUE ensemble 

We have discussed in earlier publications the possibility of computing 
topological invariants relative to Riemann surfaces using a Gaussian en¬ 
semble of iV X iV Hermitian random matrices with appropriately tuned 
external matrix sources A. The method relies on two basic ingredients 
: i) a totally explicit formula for the K-point correlation functions for 
arbitrary given source matrix A, based on the HarishChandra-Itzykson- 
Zuber integral over the unitary group [1], ii) a duality for the correlation 

K 

functions of K characteristic polynomials < fl det(AQ, — M) >^, with a 

1 

probability distribution 

Pa{M) = ( 2 . 1 ) 

Zn 

This duality exchanges the size N of the matrices with K, the number 
of points, i.e. the N x N Hermitian random matrices are replaced 
hj K X K Gaussian random matrices; the N x N source matrix A is 
exchanged with the K x K source matrix A whose eigenvalues are the 
Aq : it reads [4] 


^ n det(A„ - M)e-5tdM+A)2 

Zn J 0=1 

T N 

= f n det(a,<5„^ - 

J 


( 2 . 2 ) 


This duality is clearly well adapted to the large N limit since the r.h.s. 
is an integral over matrices whose size is independent of N. But we 
want to briefly summarize how else we have used it. 

Tuning appropriately the eigenvalues of the source matrix A, one can 
obtain in the dual version the Airy matrix model which was introduced 
by Kontsevich [8] ; it appears here as an edge singularity, reminiscent 
of the Tracy-Widom kernel [9] which governs the vicinity of the edge 
of Wigner’s semi-circle. This is done by taking for the source matrix 
A the identity matrix and considering the large N scaling regime in 
which the A^ are close to one, namely — 1) hnite, the r.h.s. of 

(2.2) becomes the Airy matrix integral introduced by Kontsevich. In 
this regime one hnds 

^ n det(A„ - M)e-5*dM+A)2 

Zn j 0=1 
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^K2^gi:|.trB3+iNtrB(A-l) 


gftrA^ 


/ 


( 2 . 3 ) 


and the r.h.s. after a rescaling B —)■ BN , (A — 1) —)■ (A — l)iV 
reduces to Kontsevich Airy integral ; 

Zkt = J (2.4) 


In the r.h.s. of (2.3) tuning the Oq, one may generate the Airy matrix 
integral, whereas the l.h.s. is still a Gaussian integral whose correlation 
functions are known explicitly. Indeed the one-point function with the 
probability weight (2.1) is given by [4] 


U{a) = (tre"^) 


1 / fr ^ ~ ^ 

a J 2i7r u — ac 


(2.5) 


in which the Oq, are the eigenvalues of the source matrix A. This formula 
is exact for any N. For the s-point function [4] the result is an integral 
over s complex variables 


U{ai---as) = (tre"i^---tre"^“) 


^CTsMv 


1 

= (^2 2-1 l^i 


S N 


n^^^^^det—^—IT n(i+ 

1 2m u, +a,- Uj 


(Ti 


Hi Ojq 


(2.6) 


Let us illustrate on the simplest example how one can use (2.3) and 
(2.6), for the one-point function. We take Aa = A for a = 1 ■ ■ ■ iF. Then 
we are dealing with 

< [det(A - M)f >A=< [det(l - iB)f >a (2.7) 

and make use of ’’replicas”, i.e. of the identity 

lim ^-^[det(A - M)]'^ = tr \ (2.8) 

K^o K dN A-M 

Since A is in the vicinity of the edge of Wigner’s semi-circle, the resol¬ 
vent has to be computed in this regime, but knowing explicitly U{a), 
this is straightforward. 

The intersection numbers of moduli space of curves are defined as co¬ 
efficients in the expansion of = tr for Zkt- The coefficients , 
the intersection numbers < >g , are dehned also as 


< r, 


ni 


■Tr,. 


>9 = 



(2.9) 
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where -ipi is called as tjj class and equal to ci(£i) with ci hrst Chern 
class and Ci is line bundle at marked point i. We have shown before 
that U{(Ti, ...as) =< tre'^i^ ■ ■ ■ tre'^®'^ > is generating function of the 
intersection numbers, since it is a Fourier transform of the density 
correlation functions, 


...as) =< J ndAie"»^»tr(5(Ai - M) >a ( 2 . 10 ) 

i 

This function provides a polynomial expansion of ai. The degree of 
total ai is equal to + |), and = 3(7 — 3 + s. This is similar 

to the evaluation of the intersection numbers by hyperbolic surfaces 
[25], by replacing the parked points by disks, whose perimeter lengths 
are /i,and the generating function of the intersection numbers is 
a polynomial of k, (i=l,...,s) and the total degree is 6(7 — 6 + 2 s [26]. 

We have shown earlier, using this strategy together with replicas, how 
to compute from there the intersection numbers of the moduli of curves 
on Riemann surfaces with one marked point [4] ; the method clearly al¬ 
lowed for more marked points. Higher multi-critical singularities, char¬ 
acterized by an integer p may also be tuned from appropriately chosen 
external sources A [ 2 , 3]. One may obtain thereby a generalized p-th 
Airy matrix model [ 6 ], taking A = diag(ai,..., Oi,...., ap_i,...., ap_i), 
with (p— 1 ) distinct eigenvalues values, each of them being times 
degenerate. The conditions[2] are 


p-i ^ 


0=1 


p-1 p-1 

= 0 {rn = 3,...,p), ^ 


Q=1 


a=l 


^ 0 ( 2 . 11 ) 


and we obtain the p-th degenerated Airy matrix model. The correlation 
functions are known in integral form 


U{ai, ...,an) 


I I -e “a 

27r 


X 


ndet( 


Ui — Mo -|- (Tj 


( 2 . 12 ) 


The one point function, which corresponds to one marked point, be¬ 
comes in the scaling limit 

[/(a) = — f (2.13) 

Na J 2Tii 

with c = AEall The Airy matrix model corresponds to the 

P Q-Q; 

p = 2 case. It is also possible to continue the model to negative values of 
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p. In particular the case p = —1 provides a generating function for the 
orbifold Euler characteristics of surfaces with n marked points, allowing 
us to recover through this method the classic results of [10, 11, 6]. 

• Lie algebras of classical groups 

The previous duality (2.2) extends to Lie algebras of the classical 
groups, such as antisymmetric real matrices for the orthogonal group 
0{2N). 


k N 

< J[[ det(AQ, • I — X) >yi=< J[[ det(a„ • I — F) >a (2-14) 

a=l n=l 

where X is 2N x 2N real antisymmetric matrix (X* = —X ) and Y 
is 2k X 2k real antisymmetric matrix ; the eigenvalues of X and Y 
are thus pure imaginary. A is also a 2N x 2N antisymmetric matrix, 
and it couples to X as an external matrix source. The matrix A is 
2k X 2k antisymmetric matrix, coupled to Y. We assume, without loss 
of generality, that A and A have the canonical form : 



/ 0 

Cli 

0 0 





— (2i 

0 

0 0 




A = 

0 

0 

0 02 

0 


(2.15) 


0 

0 

-02 0 

0 




V ■■■ 






i.e. 







A = aiv 0 ■ 

••0 

UjyV, 

V = ia2 = 

(-1 ;)■ 

(2.16) 

A is expressed also as 








A = 

Aiv 0 

• • • 0 XkV. 



( 2 . 17 ) 


The characteristic polynomial det(A ■ I — X) has the 2X roots, 
(±fAi, • • •, ±iA„). The Gaussian averages in (4.1) are defined as 


< • 

•• >A = 

1 



< • 

•• >A = 

1 


(2.18) 


in which X is a 2X x 2X real antisymmetric matrix, and Y a. 2k x 2k 
real antisymmetric matrix ; the coefficients Za and Z\ are such that 
the expectation values of one is equal to one. The derivation relies on 
a representation of the characteristic polynomials in terms of integrals 
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over Grassmann variables, as for the U{N) case, but it is more involved 

| 12 ], 

Here again the Harish Chandra formula leads to explicit formulae for 
the correlation functions. The one-point function for instance is 


U{s) 


±r±{r ( v^ + al ] 

Ns J 2m \{v + f)2 + al) v + l^ 


(2.19) 


and higher point fnnctions are also known explicitly. Therefore one 
may repeat the same tuning plus duality strategy in this case, leading 
to the desired topological nnmbers for non-orientable surfaces generated 
by these antisymmetric matrix models. 


• Superduality 


Consider 


^p,q{K • • 'h/? • • •) 


1 , nLidet(A^-M) 

nLidet(/i;3 - M) 


with 


( 2 . 20 ) 


{0{M))a = ^ [ dMC>(M)e-^*''^'+*''^^ (2.21) 

Za J 

For instance the average resolvent is given by P = Q = 1 , after taking 
derivative with respect to A and setting A = p. 


Let us recall standard dehnitions for supermatrices: let 


X = 


a a 

(3 b 


( 2 . 22 ) 


in which the matrix elements of a and b are commnting nnmbers, those 
of a and [3 anticommnting. Then the snpertrace 


strX = tra — iib 


(2.23) 


ensnres the cyclic invariance. The superdeterminant is given by 


sdetX 


det a 

det(6 — /3a~^a) 


det(a — ab ^(3) 
det b 


based on the integral 


J dOdOdxdxe''^^^ = (sdetX) ^ 


(2.24) 


(2.25) 
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where 


(2.26) 


$ = (^) $ = (x n 

The formulae are obtained either by integrating hrst the commuting 
variables, or the anticommuting variables first. We use the conventions 

W 2 = 0102 (2.27) 

and _ 

e = e (2.28) 

Finally the usual bosonic formula still holds here, namely 

str(logX) = log(sdetX). (2.29) 

We are now in position to derive the duality formula for (2.30) which 
we hrst write in integral form as 


I- N Q P 

fp,o(A„ ■■■ ••■)=/ n n n dx-jx-jeidei 

d a=l 0 = 1 13=1 




(2.30) 


or, introducing the {Q + P) x {Q + P) diagonal matrix A made of /i/j , 
(3 = 1- ■ - Q and Aq, , a = 1 ■ ■ ■ P 


N 


Fp,q{K ■■ - h/? •••)=/ n dxldxldO^pde^pie- (2.31) 

Since 


a=l 


we have 


Then 


(gtrXM)^ = 

P+Q _ 

Xha = ^ 

q:=1 


N P+Q 

tr(AX) = E E 'S’SI’S 

n=l a=l 

in which the a„ are the eigenvalues of A, 

N P+Q 

‘■■X" = E E 

a,b=l q,/3=1 


(2.32) 

(2.33) 

(2.34) 

(2.35) 
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Let us define the matrix L, {Q + P) x {Q + P) 


N 


r 


(y,l3 — 


i: KsK? = 

a=l 


Li La 

Lg 


X ■ X 6 ■ X \ 

X ■ 9 9 ■ 9 ) 


(2.36) 


This matrix Ti is Hermitian but Ts is anti-Hermitian. 

To express trX^ in terms of the matrix T some commutations are re¬ 
quired and one obtains easily 

trX^ = I] = str(r2) (2.37) 


Therefore 

Fp,Qi^a • ■ ■ h-/? • • •) 


N 


= Y\. dx^dx^d9^pd9^ 


a=l 


o- a„ ,|n$n+istrr 


(2.38) 


The SUSY Hubbard-Stratonovich transformation reads 

J dAe^^d-h^^+m = e^trr2 ^2.39) 

in which A is (P -|- Q) x (P -|- Q) and like T as far as hermiticity is 
concerned. Then 


f f ^ 

Fp,q{K dA \{ dxldxld9^9l 

a=l 

■i'S^Sp-l^trA^+strAr 


(2.40) 


One can integrate out on the x’s and 6*’s. The quadratic form in the 
exponential is 

N P+Q _ 

n=l a=l 

in which P^s = 0 for 1 < /9 < P or P^ = 1 for (P -|- 1) < /9 < (P -|- Q) 
The integration then gives 

N 

nsdet "^[(A„ - a„)(5o^ - A„^(-l)^'3] 

1 

Therefore we change Aq,^(—1)^8 —)■ Aa/s and one verihes that 

strA^ = strA^ 

8 


(2.41) 



. Then one ends up with 


Fp,q{K • ■ - /i/? • • •) = / f[sdet ^[(A^ - a„)(5c,/3 - A„^] 

(2.42) 

The above identity (2.42) relates an ordinary integral to a super matrix 
integration. In this sense it is not a full duality although it can be used 
for the large N-limit or for a super-generalization of the Kontsevich 
model. However a full superduality has been derived by Desrosiers and 
Eynard for expectation values of ratios of super-determinants [13] and 
our identity appears as a simple limiting case. 


• Arbitrary /3 

An extension of the GUE duality (2.2) to the three classical Gaussian 
ensembles GOE, GUE, GSE with respectively P = 1,2,4 has been 
derived by Desrosiers [14] , but it exchanges (3 to 4//9. However the 
lack of HarishGhandra formula for integrating over the orthogonal or 
symplectic group does not allow one to compute explicitly the /c-point 
functions and we cannot repeat the steps that we have followed for 
(3 = 2. However we have used supergroup methods to obtain the one 
and two-point functions [5, 15, 44]. 

3 GUE 


X e U{N) 


We now use the duality formula for computing the one-point function 
U (a) for some symmetric spaces. 

(i) Let us hrst consider the Hermitian case: 

From the duality formula (2.20), we obtain 




X 


det(A - X) 
det(/i — X) 


+ N(\ — fi) 


g-V(/i2_A2) 

27r^ 


c- z 




N 

dtdu [Q ( 

a=l 


Qjct — 
Qjq^ ~\~ U 


1 


u — it 

( 3 , 1 ) 


The density of state p(A) is [15] 
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^ 1 f dt j du Yj ^ dg it ^ 1 {o n\ 

N J 2ni f 2ni dg + u u — it 

By tuning the external source dg as (2.11), and taking Fourier transform of 

P(^), 

U(„) = — 1 (3.3) 


Na J 2m 

which is identical to the previous expression (2.13). 

(1) p=2 

The integral (3.3) becomes Gaussian (c=l), 

. 1 /■°° du „ 


3 r 

U{a) = -e-T2 / 

(7 J— 


—oo 271 


1 

-\/—e 12 

27ra V a 


This may be expressed as a modihed Bessel function Ki{z) 


and we have 


= Vs'- 




This explicit representation gives the intersection numbers fro Riemann sur¬ 
faces of genus g, 


< Tm >a= 7-^-T, (m = So — 2) 

® (24)s^! ^ ^ 


(2) p=3 
Then 


U{a) = - / 

^ ^ aJ 27ii 


or changing u to u = 


1 f dv 


U{a) = — i —(3.9) 
3a J 2711 

The path integral may be divided into two integrals on the real axis, above 
and below the cut : 


U{a) = Ui{a) + Un{a) 

Uj{a) = ^ 

3ai Jo 271 

Ujj{a) = 

Saz Jo 71 


(3.10) 
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Uii is the complex conjugate of Uj, and U{a) is real. The integer powers of 
a, i.e. a"', cancel. This corresponds to the spin j = 2, since a"”' 3“ = 

This cancellation means that there is no Ramond term in U{a), and only 
Neveu-Schwarz types exist. 

We have for p = 3 


U{a) = (^)[ 


Slllg. 


TT 






za 


This may also be written as an Airy function Ai{z), 


1 




(3.11) 


1 

33(73 


1 du 


1 r°° du : 






In deriving (3.11) we have made use of 


(3.12) 


poo jr j /y» Onn‘^1 ^ 

cos(t -a(t)dt=-^-[Ji/3(^^) +J-i/3(^^)] (3.13) 

The Airy function Ai{z) was used for the case of two marked points for p = 3 
in [7]. 

We obtain thus the explicit expression for the intersection numbers 

1 r(^) 


< TnJ >5 — 


(12)3^! r(V) 


(3.14) 


with n = {8g — 5 — j)/3. This condition comes from general constraint for 
the intersection numbers of s-marked points of the moduli spaces of p-spin 
curves [28], 

(p+l)(2^-2 + s) = P^nj + ^ji + s (3.15) 

i=l i=l 

The result of (3.14) agrees with (3.34) for p = 3. We have 


< 'Tip >g=l 


< 'T6,1 >9=3 


(3) P=4 


1 

< T3,2 > 9 = 2 = 0, 

1 ^ ^ _ 1 
31104’ ^ 746496’ 


(3.16) 


U{a) 





du 

27ri 


4 1 

g-o-n 


1 

4(7®A 




dx _3 _ t _ 1 „ 5 / 2 „ 1/2 

- X ‘2 

2Tii 


sm 


2 


47r(7®A 


dxx 


—X+^CF^>X 


1/2 


(3.17) 
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where the contour integral reduces to tow integrals above and below the cut, 
as for the p = 3 case. We thereby obtain for p = 4, 


1 _5/4 ^ 1 

= - n ' e 80^ > — { - 


C/(a) = -a 




n\ 






n=0 

e .0- |r(h + i^=/^r(h +Wr(b + .. 


dTT ^ 2 ' 32 

This one-point function is the generating function 

U(a) = i E < >» <^"+“4»-'r(l - 1^: 


( 3 . 18 ) 


( 3 . 19 ) 


with n = 1 ( 105 ' ~ 6 — j). Therefore the intersection numbers for p = 4 are 


< 'Tlfl >3=1 — 

< '^8,2 >3=4 = 


1 _ 3 

g. < 7 - 6,0 > 5 - 3 - 


77 

39321600’ 


19 

< rn,o >3=5- X04857600’ 


3 

20480’ 


(3.20) 


The exponent of the integrand (3.3) may be expressed as the Chebishev 
function Ti{t,x) = + Axt^ + 2x'^ with x = and we obtain a closed 

formula from the result of Appendix IV, 


U{a) = 


2^/8 


-3_c.5 

ei6o 


2sin( 


■[/- 


32 




+ /i( 

4 


32 




(3.21) 


Expanding above expression of the modihed Bessel function Iu{z) for small 
a, we have 


U{a) 


- y —=— 

8 m\n\T{n + 

1 “ 1 

8 n^=o m\n\T{n + 

1 „o3, 5 3 

—r(-)a4 H- 

87r 6407r 


T60^ ^64^ 


T60^ ^64^ 


„ / 1 6 15 

r(-)a4 + ... 


2n- i ^5m+10n-1 
2n-|- i ^5m+10n+ i 


(3.22) 


with n = l (105 — 6 — j). which agrees with the previous results [6] and 
also it agrees with the result of Liu and Xu derived by the recursion formula 
from Gelfand-Dikii equation [35]. We have obtained a closed analytic formula 
for the intersection numbers of p=4 spin curves of one marked point in (3.21) 
for arbitrary genus g by expressing it as a Bessel function. 

(4) p = 5 
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We have , 


U{a) = 

a J 27ri 


1 f du 


a J 2TTi 


dx 


-X 


_4/5^-X-fal2/5a;3/5_l ^24/5^1/5 


5(7®/® / 27ri 

By taking paths around a cut, similar to p = 3,4 cases, we have 


(3.23) 


rr/ N ^ — 


(ix^_4^5^_^|^12/5g6iri/53,3/5_^^24/5g2iri/53,l/5 


'0 27r 


1 da;^_4/5^_3,5g^l2/5g-67ri/5^3/5_^g.24/5g-27ri/5^1/5 


7o 27r 

Stt ^5^ hTT ^5^ 7207r ^5^ 

sinf 341 ^,3, 42 
5 :r(-)a^ + --- 


TT 207360 '5' 


(3.24) 


We obtain 


_ 1 _ 11 

< r.,„ J, < T,,, — 


< ^8,1 >g=4 — 


341 


25920000 
which agrees with [ 6 ] and [32]. 


, < Bo,3 >9=5— 


< '^5,4 >9=3— 0, 
161 


777600000’ 


(3.25) 


We use u = sinh 6 *, and note that T^{iu) = icosh56', 

[2 r°° 5 ll\/2 

U{a) = \ - / (76* cosh6'exp[—2a;2cosh56* H-(T®sinh 6 *] (3.26) 

V 3 JO 16 


with X = i(^) 5 (T 5 ^. By the change of 6* —)■ we have 


1 [2 1*00 


—2x7cosh6 


E 

n=0 


1 ,11^2 


n!' 16 


a‘'sinh-)^cosh- (3.27) 


This integral is evaluated by the formula. 


roo 1 

/ =- 

/o sin UTT 


r cos uO - —luiz) (3.28) 

Jo sm Z/TT 


where Iu{z) is modihed Bessel function. The genus one {g = 1) term of this 
series becomes 


U(a) ~ -\-Ki( ^ ~ -atr(l - -) + 

5 V 3 «^2y2-3i ^6^5^ 


(3.29) 


13 



which gives | for the intersection numbers of the moduli space of p = 5 spin 
curves. 

We obtain from the equation of (3.27), the intersection numbers < >g, 

with condition 6{2g — 1) = 5n + j' + 1, for p = 5, 

(5) general p 


U(a) = i / * e 

a J 2Tii 


— auP 


p{p - 1 ) 


p(p- l)(p-2)(p-3) 


X exp[-flH^^3^p-2 _ _ .. .p_3Q) 

By choosing a integral path around a cut, 


U{a) = Re{- 


e p 


pa PTT 


3!4 


We have 
U{a) = 


p(p- l)(p-2)(p-3) 4+4 2«(1-1) 1-1 
-——- a PC ^ p'x p—■■■| 


1 /. 27r 1 p -1 1+1 27r 1 

-T^(sm— )r(-) + psm— )r(l--) 

i+- p p 247r p p 


(3.31) 


npa 

(p-l)(p-3)(2p+l) ++ fa _ 3 

2 76071 ^ p’ ^ p^ 

(p - 1) (p - 5) (1 + 2p) (8p2 - 13p - 13) 5+5 . IOtt 


71433^77 


a "^p (sin- )r(l-) 

p p 


(p - l)(p - 7)(1 + 2p)(72p4 - 298p3 - 17p2 + 562p + 281) 7+1 

- a p 


9!4415 


, 147r,^, 7, 

X (sm -)r(l-) + ■■■ 

p p 


(3.32) 


The intersection numbers of p spin curves are obtained with the condition 
(p+l)( 2 p-l) =pn + j + l. 


U{cr) =J2< ^n,j >g P® p' r(l - ^^) sin — 

V p p 


(3.33) 


with m = 2 ti + 4:n{g — 1 ). 

p — 1 


< Bp >g=i— 


< Tij > 3=2 — 


24 ’ 

(p-l)(p-3)(l + 2p) r(i-^) 

i+> 


p5!423 


r(i - 
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< TnJ >g=3 


^ ^ g=A 

r(i-p 


{P - 5)(p - 1)(1 + 2p)(8/ - I3p - 13) r(l - |) 

^2714332 Y {1 - i±i) 

{P - 1)(P - 7)(1 + 2p){72p^ - 298p3 - I7p2 + 562p + 281) 

p39!44l5 


(3.34) 


This result is same as [ 6 ], where the integral is restricted to a path from 0 to 
cx) without sin ^ factor in U{a). 

(6) p = -1 

This expression for arbitrary p in (3.34) allows the analytic continuation 
to the negative values of p. In the case p = —1, it correspond to Euler char¬ 
acteristics x(Mg^i) = C(1 — 2g) [ 6 ]. For p = —1, the power of becomes 
zero, and a dependence disappears. Therefore, we need the introduction of 
c, which is taken as N to specify the genus g. 


rrr ^ I f du -Nlog^ 

Na J 2tii 
_ 1 f du u — ^ 


N J 2TTi n + I 

POO 1 R 

•^0 1 - e * 


2nN‘^^ 


(3.35) 


We have used 7 ^ 4 - = e and a following expansion. 


1 - e-* t 2 


1 1 OD f2n—l 


n=l 


(3.36) 


This gives Euler characteristics (intersection number < t >g for p = —1), 

X(M,,0 =< r >,= = C(1 - 2g) (3.37) 

where ( is the Riemann zeta function and is a Bernoulli number {Bi = |, 
B 2 = ^, R 3 = R 4 = ^). When p is negative, we have to specify the 
meaning of spin j. This index p is related to the level k of the Lie group 
su{2)k/u{l). This was studied by Witten [36] as a chiral ring (Landau- 
Ginzburg theory) of primary helds and their gravitational descendants with 


p = k + 2 


(3.38) 


The present case corresponds to the singularity theory of Ap_i. When p 
is negative, we have a non-compact Lorenzian group sl{2, R)k/u{l), whose 
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discrete spectrum is known to correspond to two series and D- [37], The 
analytic continuation of p —)■ —p corresponds to and the spin j takes 
negative value. For instance, for p = —1 , the Euler characteristics x{^g,i) 
is dehned by the top Chern class only and n , which is the power of the hrst 
Chern class, should be zero. Then we have only < tq-i >g in which j = —1. 

< Vi ><;= 4 / C't(z/)[ci(£)]" (3.39) 

Thus < Tn-i >g is not surprising since the discrete spectrum with negative 
spin exists for SL{2, R). 

(7) p = -2 

As noticed in [38], we have two expansions, weak conpling and strong 
conpling for p = —2, which correspond to the Gross-Witten model for the 
nnitary group [27, 39]. There is a phase transition between these two phases. 
The weak conpling corresponds to small a and strong coupling corresponds 
to large a. Therefore, the spin values j takes negative valnes of for weak 
coupling, and positive values for the strong conpling phase. The expansion 
of u(a) is expressed by putting n = 0 as 

l_|_j 

u(a) = ^ cija~ (3.40) 


More details of the discrete spectrnm of SL{2, R)k are presented in appendix. 
Before closing this section on GUE, we write the one point function U (a) as 
an angular integral, which is useful for the strong coupling expansion. 

In the expression of U (a) one pnts sin 6 = l/vTT”?? and cos 6 = u/\/l + 
in (3.3). Then with a = it, and u = |n. 



dv 


exp[ 


, t 


,p+i 


P + r 2 ' 


{(n + i 


\p+i 


With V = it becomes 

sint^ ’ 





1 I 2ic i^^, sin(p+l)g 

27r (sin Oy p + 1 2 (sin 6 ')p+i 


(3.41) 


(3.42) 


Note that the denominator of the exponent (sin0)^+^ becomes a numera¬ 
tor when p-|-l is negative, and it provides a large a expansion (strong coupling 
expansion for large t) corresponding to a discrete spectrum of SL{2, R)k. This 
large t = —ia expansion becomes, for instance for p = —2, 


U(a) 


yiD-yny- 

^ m ^ 



7 ,c,. 

—(-) + 


(3.43) 
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with 


(3.44) 


(2m. - 1)! 1 

"■ m! nE'(-P) 

D is a divergent term, which should be regularized. The above expression 
matches exactly a strong coupling expansion for the unitary (gauge) group, 
for a single trace result with iV = 0 [38]. The unitary matrix model is 

Z = J df/e*’^(uct+utc) ^3 45^ 


where U is a N x N unitary matrix, UW = 1. C* is an external complex 
matrix. The strong expansion is an expansion in powers of tr(ClC)™. The 
coefficient of tr(ClC)™, Cm is equal to 


’^2- JV2-1’ (iV2-l)(iV2-4)’ 

30 

“(iV2-l)(iV2-4)(iV2-9)’'" 


(3.46) 


For obtaining the N dependence, we need the insertion of a logarithmic 
term in (3.3) as [38], 



u + 1^ 


(3.47) 


4 Classical Lie algebras 

X g 0{2N) 

When the random matrix X varies over a classical Lie algebra, with 
Gaussian distributetion, the n-point correlation function in an external source 
is obtained again exactly , after use of the Harish Chandra formula[43]. We 
have discussed in earlier work such models with external source[12, 44]. 

Consider the Lie algebra of 0{2N), namely real antisymmetric matrices. 
Since, the Harish Chandra formula holds for this Lie algebra, we can obtain 
explicit expressions for the n-point correlation functions. Again one can 
derive a duality identity. In the present case, instead of the duality formula 
involving a supermatrix Q, it is convenient to use 

k N 

< n det(A« • I - X) >A=< n det(a„ ■ I - X) >a (4.1) 

a.=l n=l 

where X is a 2N x 2N real antisymmetric matrix (X* = —X ) and Y is 
2k X 2k real antisymmetric matrix ; the eigenvalues of X and Y are thus 
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pure imaginary. The matrix source A is also a 2N x 2N antisymmetric 
matrix. The matrix A is q 2/c x 2/c antisymmetric matrix, coupled to Y. We 
assume, without loss of generality, that A and A take the canonical form : 


A = aiv 0 • • • 

eUATU, V = ia2= Q 1 . 

(4.2) 

A is expressed also as 

A 

= Xiv © • • • © XkV. 

(4.3) 

The dehnition of the averages are 


{0(X))^ = 0 

J dX 0{X) exp (^trX2 + trAX) 

(4.4) 

( 0 (K))a = 0 

0 A 

J dY 0{Y) exp (WrY^ + trAY) 

(4.5) 


By an appropriate tuning of the a„’s, and a corresponding rescaling of Y 
and A, one may generate similarly higher models of type p with the conditions 

( 2 . 11 ), 

Z = j (4.6) 

where p is an odd integer. 

The HarishChandra integral for the integral over g G SO{2N) group, and 
given real antisymmetric matrices Y and A, reads 


N 


'Z (dettc)exp [2 X! 'w{yj)Xj] 

>soi2N) ^ n iy]-ym-\l) 


(4.7) 


l<j<k<N 


where C = {2N — 1)! n|=r^(2j — 1)!, and w are elements of the Weyl group, 
which consists here of permutations followed by reflections {ui —)■ ±|/j ; i = 
1, • • •, N) with an even number of sign changes. 

For the one point function, we obtain when X is a 2N x 2N real anti¬ 
symmetric random matrix, from the above formula. 


U{a) 


2N 
2N 


< tre®"^ 


Na 


N (n zY 

E n( " 

a=l 77 ^ 0 : 

du .{u + |) 
27ri 


a 2 -a 2 
2 ^ 




U 


U 


4 - 1 - (a ^ —a) 




(4.8) 
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where the contour encircles the poles u = a^. Or, shifting u ^ u — j, 


U{a) = 


dv 


{u - 1)2 - a? 


u — 


Na J 2m “ {u + |)2 - aj 


4 g^n 


U 


Tuning the external source to obtain the p-th degeneracy, one hnds 


1 




du 


-^[o+ff+'-o-fr+h 


Na J 2m 


a 


il- — 
^ Au 


(4.9) 


(4.10) 


(l)p = 3 

There are two terms in (4.10) ; the hrst term U{a)^^ is exactly one-half 
of 17(f) for the GUE (orientable Riemann surfaces). The second term is a 
new term, and we denote it as the non-orientable part U((j)^‘^, since it is 
related to non-orientable surfaces with half-integer genus : 


U{a 


,OR 




2-33(1)3 


Ai(- 


■4.31/34^"' 


(4.11) 


For the non-orientable surfaces, from the condition, 
(p + l)(2p- 1) =pn + j + l 


we hnd that the genus g is always a half-integer [g = 
has a series expansion in powers of a^~^ 

U{af^ ■- 


13 5 
2 ’ 2 ’ 2 ’ 


(4.12) 
), and U((t)^‘^ 


1 1 

du 1 

AJ 

2m u 

1 

dx 1 

12 ] 

2m X 

Re{ 

^ / 
-1 • / 


i+j 

p . For p = 3, we have 


e 

■CO 


(4.13) 


This function may be expanded as 

1 


U{af^ = Re{- f 

^ ^ ^ 12m Jo 


12m 
1 1 a, 
7r^^2' 


00 1 

dx—e~^ 

X 


OD -1 I 

n=0 ^ ^ 


271, 

3(sm — )r(-; 


,a 


3 ' ^3' 

1 1 /CT 16 . 47r 2 1 1 

^^*^2^ ' ~ ^3! ■ 12 -43'2 

1, ^ r^)f(sing^)r(l)---- 

7r4! ■ 12-444^ ^ 3 ^ ^3^ 


-)«(sin27r) 


(4.14) 
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Using Airy functions, the p = 3 case is expressed as 


U{a) = 


2.31/3(1)4/3 


1 

Ai{x) — - dx'Ai{x') (4-15) 

4 Jo 


with I = 

The Airy function Ai{z) and the integral of Airy function may be ex¬ 
panded as 

(4.16) 

^onlT{n + ^) 3 3^i'^ ^ n\T[n + ^) 3 


At{t)dt 32/3r(|)^ 34 / 3 .2r(|)^ ^ 36-32/31(1)^ +■■■ 

Inserting these expansions, we have for p = 3, = ^(1)^(1)- 

Uia) = _ - __ - _(-)20/3^... 

/ ' O/i'P/'l^ / o / mo eiATf2\\ o' ' 


24r(4)^2 


108-64r(|)^2' 


^ ^48r(f)^2^ ^ 384r(l)^2 


5_/^A6/3_!1^^_|'^'|32/3 


864-441(1)'2 


+ ■ ■ -I 


(4.18) 


U(a) = <T,,„>,.,r(l-i)(|)‘+i+<T«,,>,.3r(l-2)38(|)«+i + ... 

+ [ < t2,l >9-3/2 r(l - ^)3^( —< Ts.o >9-5/2 T)! — g )8'^(^) 

+ < Tio,. > 9 . 9/2 r(l - /)3*(|)88/8 + . . .] (4.19) 

We have for p = 3, 

Uiapo ^ J_,2r(i _ 2, ^ _ 1, 

+^i/»r(l-H)4.... (4.20) 

We have obtained for p = 3 the explicit intersection numbers for non- 
orientable surfaces with one marked point. The intersection number < 
^ 2,1 >g= 3/2 corresponds to a cross-capped torus. For p = 1/2 wa are dealing 
with the topology of the projective plane but for this case, the intersection 
numbers < Tq i >g=i /2 are present only beyond the two marked points level 
[12]. We have 


<C Ti 0 1— -, To 1 — -, 

9-i 24’ 2,i g-- gg^, 


(4.21) 
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(2) general p 

Using the binomial expansion, one finds {y = 

U{a) = - ^ H-] 

^ ^ AypN J ^6 ^ 

x[l + yrp] (4.22) 

This is again the sum of two contributions, orientable (OR) and non- 
orientable (NO). The odd powers in y correspond to the orientable contri¬ 
bution, which is the same as for the unitary case ; the even powers in y 
correspond to the non-orientable case : 

U{a) = U{af^ + U{a)^^ (4.23) 

U{a)^^ is same as GUE but the normalization of a is replaced by a 12. 

The first term in the above series expansion is divergent, and it should 
be regularized. Except for this divergent term, we give the series expansion 
up to order y^ (we have neglected the phase factor sin(^^), 

uwf° = f^(p-i)r(i--) 

24 p 

+ | 7 (p-- 5p+i)r(i --) 

o! p 

+ ^{p - 1)(P - 3)(4p3 - 23p2 - 2 p - 6)r(l - ^) 

+ ^^j^|^(p-l)(9p6-l21p5 + 435p^_317p3 

-167p2 - 471p - 43)r(l - -) + Oiy^^) (4.24) 

p 

From this genus expansion, one obtains the intersection numbers of p-spin 
curves for non-orientable surfaces. 

(3) P = -1 

We now perform the limit, p —)■ — 1, which is related to the virtual Euler 
characteristics. When we put p = — 1 in (4.24), the T function term becomes 
an integer for p = —1, and this agrees with the intersection number of < 
U,o >g, which gives a factor r(l ~ “) = r(2) = 1 for the spin zero. We 
obtain 

U{a)^^ = - —(2p)2-—(2p)^-—(2p)6-—(2p)S + ---(4.25) 

V ; 24^ 240^ 504^ 480^ ^ ^ 
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This series agrees precisely with the series expansion 


OO 1 1 

= -1: (4.26) 

g=l 

where Bg is a Bernoulli number, a positive rational number. i?i = |,i ?2 = 
^,i ?3 = ^,i ?4 = The coefficient of is the same as for the virtual 

Euler characteristics of the moduli space of real algebraic curves for genus g 
and one marked point, which was derived from the Penner model of the real 
symmetric matrix by Goulden et ah [42], ( We use for the half genuses in 
the list , 1, |, 2,... for a projective plane, Klein bottle, cross-capped torus, 

doubly cross-capped torus with the notation g = 1, g = 2, g = 3,g = 
4,....,respectively [41], and this is a reason for the appearance of the (2?/)^^ 
factor in (4.25)). 

Since we derived this from the antisymmetric 0{2N) Lie algebra, the 
coincidence between 0{2N) he algebra and GOE for the virtual Euler char¬ 
acteristics seems remarkable. 

= (4.27) 


This result may be obtained analytically to all orders. We now derive this 
result from the integral form (4.10) replacing c by N. With p = —1, it 
becomes 


U{(t)= ,1, f duC 

4Na J + (j' u 

(4.28) 

With the change of variable u —)■ au, 



(4.29) 

We divide it into two parts, and , 



(4.30) 


(4.31) 


We use the same change of variables as for the unitary case [4], 


?y — 1 1 p ^ p ^ 

- -= e-y, u = , du = -2-—- -dy (4.32) 

u + l ’ i-e-y' {i-e-yf ^ ^ ^ 
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(4.33) 





, e-y ,l-e-y. 

(1 - e-y)‘^^l + e-y' 
1 1 
^1-e-y l + e-y' 


(4.34) 


It is interesting to note that both Boson and Fermion distributions enter 
in the above integrand (4.34). 

If we use the expansions, 


1 11 °° 


n-1 _2n-l 


^ _ 1 , ^ (-1)^-H2^^ - 1) o ..2n-l 


1 + e-y 


(4.35) 


then they become 


^ ’ 2N J ^ 2niV2« 

TTI^\N 0 _ ^ f 1 T\n-lBn 1 


^ dye-^^+iEi-^r-^^ 


2n iV 2 ^+i 


4^/ ^ 2 n ^N^n+i 


= -T^r dy- 


^y I ^_(2-22-)5„ 1 


7E(-1) 


y ^n=l 


2n iV 2 ^+i 


(4.36) 


We now get from the above equation (replacing n by 5 ^), 

= -icd - 2 ,,) = 


(4.37) 


For s marked point, the result obtained from the real symmetric matrix 
Fenner model [42] is 


X (A4,)-(-l)--- B, 


(4.38) 


This result can be obtained by applying equation (4.37) [7]. In this 0(2N) 
model, we have the following condition, the same as for Riemann surfaces 
with spin j and s-marked points 


(p+ l)(2^-2 + s) =pE^i + Eji + 


(4.39) 


2 = 1 2=1 
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However, we have to assign the genus g also to half integers to represent 
non-orientable surfaces [12], 

X € 0(2iV + 1) 

For SO{2N + 1) Lie algebra, the matrix X is 

X = hiv © h2V © • • ■ Hnv © 0 (4.40) 

The measure is V{HY, 

N 

v(H)= n (4-41) 


The Harish Chandra formula is 


N 

Y, (dettc)exp(2 Y w{aj)bj) 

- — - Jj - (4.42) 

n (oj - - bl) n ajbi 

i=i 

with Coin) = nf=i(2j - 1)! U-= 2 NJ'- Comparing with the 0{2N) case, this 
formula differs from (4.7) by the presence of the term JJajbj in the denomi¬ 
nator. For the one point function, we have 


I = 


ISO(2N+l) 


U{a) 


1 

X 


E 


a=l 



N 

i=l 


n(A^ 

n(a? 


E ^ ai\i 

a])l\ak 


(4.43) 


This sum of integrals may be written as a contour integral, which collects 
poles at u = of, 


U(a) = I 


du {^/u + aY - Oj 

O/TT/l AA ni _ 


a 


{n=a2} 27ri El 


TT ('^ Clj V + a 

27J'* — a] a + 2v 


2 
a 

1 r dv 

a f 2ni (n - |)2 


{^/u + ay 


— (1 H—^)e' 
u Ju 




j=i '' 

^ (v+^f 


n 




a 




(4.44) 


By the tuning to the p-th degeneracy, we obtain 

[/(a) = ll (4.45) 

^ ’ aJ 2^ ^ 2u’ ^ ’ 

This takes the same form as for the 0{2N) case. 
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(4.46) 


X € Sp{N) 

The Haar measure of Sp{N) is A(A)^, with 

A(A)=no?-\")n^t 

i<j k 

The Harish Chandra formula for Sp{N) reads [44] 


I = 

= C 


Jg ^ A(a)A(6) 


det [2sinh(2aj6j) 


- a])m -h])Y{Mk) 

For the one point function, we have 


(4.47) 


][ ^ poo ^ 


n ( Af - A2 ) n A , 


Uicr) = IlA 


l<i<j<N 


a=l ' 


' 2=1 


du ^ {^/u + a)^ - 


n {al - a ]) n ak 


l<k<N Y, A 2 +aA „+2 x; Ai 


l<A:<Ar 


u — (i/u + ay — u 




a‘^-\-2cFy/u 


dv 


N 


(v + ay v + a 9 


^ J gL U ~m — r 

a f 2m — aj a + 2v 


.o-^-\-2av 


dv (v + 1)^ - 




a 


a J 2m M (v — ^y — a 




(4.48) 


where we have shifted u — )■ u — | and — )■ a^/2. This expression becomes 
the same as for the 0(2A) case, when we put u ^ 2u up to a factor 2. Note 
that we do not need to consider the expansion | as in the 0(2A) case. The 
hrst term of the expression is same as for CUE. By the tuning to the p-th 
case, we have 

U(a) = i / Ae-riT«”+!)'«-(“-!)'«)(i + f.) (4.49) 

^ ’ aJ 2m ^ 2u ^ ’ 

We write these two terms as U{a) = U{a)‘^^+U{a)^‘^. It is then obvious that 
we obtain the same intersection numbers and virtual Euler characteristics as 
in the 0(2A) case. 


5 Open intersection numbers 


Kontsevich-Penner model 
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The Airy matrix model with an external sonrce, the Kontsevich model 
for p = 2, gives the intersection nnmbers for closed Riemann snrfaces, which 
satisfy a KdV hierarchy. These closed intersection nnmbers are obtained 
from (3.33) and (3.34) for one marked point. They are known for genus g 
and one marked point in a simple closed form, 

_ 1 

When the Riemann surface has boundaries, open intersection numbers ap¬ 
pear, which differ from that of the Kontsevich model. We have studied earlier 
the effect of an additional logarithmic potential in the Kontsevich model, the 
so called Kontsevich-Penner model [19]. In our work this model came from 
a two matrix model, which originated itself from a time-dependent matrix 
model. The eigenvalues of the two matrices correspond for one to the edge of 
the distribution and for the other one to the bulk. Then we can use the du¬ 
ality identity for the two characteristic polynomials of the two matrices with 
external sources, and thereby recover the Kontsevich-Penner model. There¬ 
fore the presence in that model of the term i of (detM)^ = exp[/ctrlogM] 
coreesponds to the addition of a boundary (an open disc) in the random 
surfaces described by the Kontsevich model. Recently the open intersec¬ 
tion numbers have been analyzed in [16, 17, 18]. The generating matrix 
model for those open intersection numbers are given by a Kontsevich-Penner 
model [21, 22]. This Kontsevich-Penner model has different Virasoro equa¬ 
tions and different intersection numbers, which depend upon an additional 
parameter k which corresponds to the logarithmic term 


(5.1) 


^ _ J ^j^g^trM^+trMA+ktrlogM 


For the open intersection numbers, considered by [16], k takes the value k = I 
[22]. The addition of the logarithmic potential yields new Virasoro equations 
and new intersection numbers which related to the boundary insertions. The 
intersection numbers for the model (5.2) have been computed in [19], 

1 -f 12k'^ 

< Ti >g=i= ———, < Tori >3=1= k, ... (5.3) 

The appearance of a half-integer index exhibits the non-orientable nature. 
The non-vanishing < ■ ■ ■ Tn^ >g are restricted by the condition 


3{2g — 2-|-s) = 2y^ni + s 


When the parameter k vanshes, the intersection numbers reduce to the usual 
Kontsevich result, which satishes a KdV hierarchy. When the cubic Airy 
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matrix part is absent, and only the logarithmic potential is present (Penner 
model), as we have seen in the p = —1 case in section 3, the model gives 
the Euler characterstics [6]. When A; = 1, it reduces to open intersection 
numbers. The meaning of the parameter k is found in the two matrix model 
[6, 19], 

We now consider the /c-dependence with one marked point. The one point 
intersection numbers of the Kontsevich-Penner model (5.2) are obtained from 
U{a) [19] 


U((j) = - <p -e ■ 2 

^ ^ aJ 27 ri 




--£-0-3 
= —e 12 

a 


du 

2m 


^—aTu‘^+klog{u-lr^(T) — klog{u—^a) 


with 


_ 1 _ 1 

'"“A’ ‘"“W 


(5.5) 


(5.6) 


This U{a) correctly reduces to the intersection numbers of the Kontsevich 
model with one marked point when k = 0, 




vr (—cY 




t3g-2 


(5.7) 


Including the factor (p=2), the intersection number reduces to 


< '^33-2 >g- 


m^gi 

For dealing with higher /c’s, we expand (5.5), after rescaling of u, 


(5.8) 


U{a) 


1 CO- 3 r dtl c 2 [- ,,2 3 2 9 2 

—e 12 f—e 4 1 +/c(_a 2 + _a 2 + —a 
t 2 J 2m u 6u-^ ou° 


2a 

2 4 

kHA:a^ + 




3u^ 






) 

(5.9) 


The coefficients of the successive orders in k tmay be computed from 


k 


3 

(J2 


du _£„ 2 , U + 

a3/2 




3/2 

2 


) 


(5.10) 
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with 


erf(a;) = —j=^ 




dt 


TT Jo 


(5.11) 


where the integral is computed as the discontinuity across the cut between -1 
to 1 in the w-plane. This integral becomes a contour integral around u = 0 
by expanding the logarithm in powers of - as 


k C<T^ 
rC 12 


a 


3/2 


du X—^ 

-e 4 “ \ 

2711 ^ 


^§(2n+l) 


{2n + 


= ke 12 


ca'^ 


„ 3 , n'.(2n + 1 ) 


(5.12) 


For the odd powers of k, the integration over u is the same contour integral 
around u = 0 [19]. 

For the even powers of k, we use the following integrals, 


r due-^^'^ = (- 1)11 2 :^ 

1-00 M2n ^ ^ ( 2 n-l)!! 


(5.13) 


which may be obtained by integration over a. Putting a = | we obtain 

r c ,,2 1 


e 4 


:du = 




cn 


(5.14) 


Thus we obtain, up to terms of order k^a'^, 


U{a) = /—e-i“'[l + 2A;V 

~ J 2ti 


2(7 2 






3 

(T 2 


(5.15) 


Expanding the factor e 12 *^ , we obtain the intersection number < ti > as 


< Ti >3=1- —(1 + 12/c^) 

For < 74 > 3=2 and < 77 >3=3 , we obtain with (5.13) 

1 


(5.16) 


< t -4 > 3 = 2 = 


< 77 > 3 = 3 = 


1152 

1 


(1 + 56r + 16r), 


-(25 + 5508A;^ + 3120A;^ + 192A;^) (5.17) 


2073600' 

The intersection numbers for fractional genus, < 75 >„= 3 , < 711 >„= 5 , 

2^2 2 ^2 

are expressed as polynomials with odd powers of k and they are given by 


the residues for the terms of order a®,... in (5.9). 


< 75 >„_3= —(k + k^ 
2 9-2 12^ 


(5.18) 



In [19], there is a misprint for this term of order which had been eval¬ 
uated from Virasoro equations. The above results agree with the Virasoro 
equations, which will be discussed below. In general, the intersection num¬ 
bers with one marked point < T 3 g _2 >g are easily computed to all orders by 
using the formulae (5.12) and (5.13). 

We have used the condition corresponding to p = 2 and one point, s = 1, 
(p -|- l)(2p — 1) = 2n -|- 1 for < >g. This condition 3(2p — 1) = 2n -|- 1 

implies n = 3g — 2 and if n is a half integer, then the genus g is also a half 
integer. Those half integer g appear for non-orientable surfaces, as discussed 
earlier with the random surfaces generated by antisymmetric matrices [12] in 
section 5, and it corresponds to the topology of non orientable surfaces such 
as the projective plane (p = |), the Klein bottle (p = 1), the cross-capped 
torus (p = |), etc. [41]. 

The string equation for the Kontsevich-Tenner model has been derived in 
[19, 22]; it reads 


dF 

dto 


(^ + 7^)tn+l 
n=0,l,2,... ^ 


dF 

dtn 


+ E 


n=i ^ 
2 ’ 2 ’ 


. K OF 


1 2 k 

4 0 2 2 


(5.19) 


The free energy F is divided into close and open parts, F^ and F°. Then, 
we have 


dF'^ 


dF° 

dto 



E + 9)^1 

n=0,l,2,.. ^ 


dF 

dtn 


k 

2 


tl. + 



in + -)tn+l 


dF 

dtn 


(5.20) 


The Virasoro equations for open intersection theory for genus zero has been 
discussed in [16, 17]. The open intersection numbers are dehned analogously 
to the closed case as 

< Vn^ ■ ■ ■ Tn,a^ >g= f_ • • • C" (5-21) 

Jm t 

g,k,s 

The string equation becomes for the open free energy is 


dio 


E ^*+1 


i=0 


dF° 

dii 


+ s 


(5.22) 


which is consistent with (5.20) (s is proportional to k, and the difference is 
due to a different normalization of ). The string equation implies 


< 


roll 


^ \ ' 


:=E< 






a >, 


(5.23) 
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We will consider the case of two marked points and derive this string equation 
in the next section. 


open p-th spin cnves 

The open intersection numbers for the p-th spin curves with boundaries 
are also given by the addition of a logarithmic potential to (3.3) 


^ ^ -^[(«+fy+l-(n-f y+l]+fclog(n+f)-fclog(n-f) 


U{a) = — f —e p+i 
(j J Stt 

Expanding the exponent, 

U{a) = -<f 

a J 2tii 

_ p(p-i)(p-y)(p-3)^5.„,-4 


(5.24) 


3!4 


+ 


a 


5!42 

■) + IkH-y + W-* + 




x[l + A;(-a+- , 

^ M 12^3 80^5 

+■■■)++■■■)+ o[e)] 

3! m3 4! 




(5.25) 


By choosing an integration path around the cut, with x = the above 
equation becomes 


U{a) = 


i+- 
pa PTT 


dxx 


- —1 —X 

ve 


^ r p(p-l) 2+1 1-1 p(p- l)(p-2)(p-3) 4+4 4_4 
xexp[-; 7 ;-;- a PX P -^777^- (J PX P + 


3!4 


5!42 


r-. 7 / 1+i -i 1 1+1 -1 1 1 

[1 + k(a ^px p H- a^^px p H-jo 

12 5*2 


’X p + • ■ ■) 


^2 

2 

1 


1-1 1 44_4 _4 1 PI 3 _3 

— {a"^px pH— a^px p + • • ■) H— -k^{a^^px p H-■ ■ ■) 
2 6 3! 


+ ^k^{a^"^px pH-) + 

The integration over x gives 
1 k\ 1 


(5.26) 


U{a) = -(^—^ + —)-a^+pr(l - -) - i—k + —A;3)-a^+lr(l - -) 
^ ^ ^ 24 2 n P 24 12 n ^ p^ 


144^ 


(p-l)(p-3)(l + 2p) 
40 


+ (3p + 1)^ + 2A;^]-a^+ir(l - -) + 

71 p 


(5.27) 


This expansionc provides the following open intersection numbers, 
p — 1 H- 12A;2 


< B,o >g=i— 


24 
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< r2^ >g=i= ^{pk + 2k^) {p^2), < r| >^^1= ^(/c +/c^) {p = 2) 

1 ip-l)(p-3)(l + 2p) 2 o;.4i 


< TnJ >g=2 — 


p{i2y 


40 


- (3p + i)r- 2 r 


r(i-i^) 

(5.28) 


where the condition [p + l)( 25 f — 1) = pn + j + 1 determines n and j. For 
p=2, g = 2, we have from the above expression, 

1 

< ^4,0 > 9=2 - 


-[1 + 56r + ler 


(5.29) 


(24)22!' 

which agrees with the result of (5.17). The higher order open intersection 
numbers of p-th spin curves and one marked point are easily evaluated from 
the expansion of (5.26). The intersection numbers are related to U{(7) as in 
[6] 

= >g r(l - (5.30) 

For p = 2, (Kontsevich-Fenner model), we have 


U{a) = —12^^^ lim [ 
2tT(T^ P^2Jo 




dxe ^xp 


1 , 1 i+i ^ 
1 + px p 

1 1 i+- -- 

i — PX p 


(5.31) 


Additional computations of open intersection numbers for p-spin curves are 
listed in an Appendix. In this Appendix, we derive also the string equation 
for p spin curves in the presence of a logarithmic potential. 


open 0{2N) model 


We now consider the non orientable intersection numbers provided by 
the 0{2N) model. It is natural to investigate the relation between the non- 
orientable intersection numbers given by the 0{2N) model and the open 
intersection numbers which we have just discussed. The open intersection 
numbers for the 0{2N) case with a logarithmic potential is also interesting 
since the model deviates from KdV and KP hierarchies. 

For the 0(2N) case with a logarithmic potential, U{a) for the p-th higher 
Airy singularity becomes (4.10), 




1 a 


2TTi 




(5.32) 


Since it resembles to the unitary case, with the replacement a by 2s, the 
expansion (5.26) can be used. 


U{s) = 


pS PTl 


1 1+1 _i, 

dxxp e (1 — -s PX p^ 
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p{p-l) 2+2 1-2 p{p - l){p - 2){p - ?>) 4+4 4_4 

xexp[ -^ 77 -; - S J>X P - 7777 - S PX P+---| 


3!4 

i+i -2 


5!42 


lc:55 

[1 + k{s'^~''px~p + —s'^~'~px~p + -— ^s^~^px~p + • • •) 

LZi 0 * ^ 

k'^ , 2+^ 1 4+^ \ 1 ,S/ 3+^ 

+ —(s PX p + -s PX p +■■■) + -;:Tk (s PX p + 


2 ' 6 
-^k'^(s^~^px~p + •■■) + •••1 


3! 


(5.33) 


1 1-I-- —-\ 

The term (—^s px p) gives an additional contribution to the open inter- 
section numbers characterized by a parameter k as discussed in (5.28). This 
contribution reads 

U{s) = Uo{s) + AU{s) 

Ar7'/'^ 1|.^ V I P ^ 2+--p^.| 2 'j 

Ac/M = ,r(i.r(i--) 


(5.34) 


where Uo{s) is the same as U{a) in (5.26). Thus the open intersection num¬ 
bers for the 0{2N) case ( 0{2N) p-th Airy matrix model with a logarithmic 
potential), together with Uo{s), are given by 

p - 1 + 12 A; + 12 A ;2 

< 'Tip >3=1— - 


24 


< r2,i >_3 = 


(p — 1) -I- 2pk + Qk"^ -I- Ak^ 


, < r| >„=3 = 


1 -|- 4A; -|- Qk^ + Ak^ 
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< >»=2= + + 2k - (1 + 3p)k--‘ + 41’ - 21’] 


3 2 48 ’ 2^2 


p( 12)2 


40 


X ■ 


r(i - 


r(i 


i+p 


(5.35) 


6 Multiple marked points and Virasoro equa¬ 
tions 


string equation 

The Virasoro equations have been investigated for the Kontsevich-Penner 
model [19, 22]. The hrst Virasoro equation, or string equation, reads (5.19) 
and (5.23). 

s s s 

< ^0 n >9= T < n >9 (6.1) 

i=l j=l i^j 
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Since the intersection numbers for s-marked points are known explicitly from 
the integral formula ioiU{ai,as), it is interesting to derive the above string 
equation and the other Virasoro equations for the Kontsevich-Penner model 
from our formulation of the s-point correlation function U{ai, ...ag)- 

2 marked points : 

The two marked points correlation function U{ai,a 2 ) is [6] 


U{ai, (T 2 ) = e 12 




dUidU 2 r 2 2 , ,1 /“l + 

——^ expj-aiMi - a 2 U 2 + fclog(-f— 

(27r*)"' ui - 20-1 


+klog{ 


U2 


1^2 


U2 


1^2 


)} 


(ui - U 2 + i(ai + a 2 ))(u 2 - Ml + \{ai + (T 2 )) 


( 6 . 2 ) 


Writing the denominators as principal parts integrals as in [45] 


1 

Ml - M 2 + K^-i + (^ 2 ) 


—i 



-“ 2+5 (o-i+ 0 - 2 )) 


(6.3) 


we obtain 


U{ai,a 2 ) = e I f dad/Sexp{-aiul - a 2 u\ 

J [imy JO 

+k\og C^ + /dog( ^^ 1^^ ) + *(a - - M2) 

Ml - 20-1 M2 - 2<^2 

+*( 0 ; +/3)-((Ti + (T 2 )} (6-4) 


Replacing Ui —)■ and a i/Tidga and jd —)■ ^/a^a^ld, we obtain 


+ 

+ 


U{ai, 0 - 2 ) = ^/aT^e 12 ^ dad/3exp{-ul - uj 


{27ri) 

i{a - /3){^/^Ul - ^/^U2) + i-{a + (3){^/a^{ai + 0 - 2 ) 


1 3 / 2 . 

Ml - ^O-l 


1 3/2- 

M 2 - 2*^2 


(6.5) 


Keeping the term of order in the pre factor and putting the a 2 to 0 in 
the exponent, we obtain 


U{ai,a2) 


du2 

27ri 







exp[—(TiM^ + k\og{ 


Ml - iai ^ 

( 6 . 6 ) 
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Using the argument of ( 5 . 13 ) for the integration over U2, we have for the one 
marked point open intersection number < r„_i >g for the Kontsevich-Penner 
model, 

< roTn >g = < Tn-l >g ( 6 . 7 ) 

which gives the results of the string equation (6.1) ; for instance, 


< TnTi >„_!=< r_i >„_!= k 

2^2 2^2 

1 _L 1 9^2 

< Tors ><;=1 = < Ti >3=1= —-- 


( 6 . 8 ) 


In Appendix, we show that the string equation holds also for the p spin 
curves. 


3 marked points: 

The three point correlation function f/((Ti, <72, <73) is given by 


U((7i,( 72,(73) = e 12 


= I J([ n 


^ 7/ • -U -rr- 


1 


i=l 2*^* 


[Ui — U 2 + |((7i + (T2)][u2 — M 3 + |((72 + (73)][m3 — Mi + ^(<73 + (7i) 


( 6 . 9 ) 


The three denominators are replaced by integrals ofver a,/9,7 as in ( 6 . 3 ). 
Changing variables m* —)■ a —)■ fd —)■ 1/0707/9, 7 —)■ 1/07077, the 

above expression becomes 


U(o-i,o-2,o-3) = e 12 


= p“i^(°'i+'^ 2 +y 




^ M* + , 
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2=1 ^2 


1 3/2- 
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-('^S+O'l)) 


( 6 . 10 ) 


Keeping 1/07, and putting the other 0^2 to zero, this expression becomes 
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27ri M 2 
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27ri M 2 
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{ 2 mY M3 - Ml + |(o-i + 0-3) 
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Ul — Uz + I(<71 + (73) / '»l + CTi ^fc^ M3 + q'3 ^fc 
ui - + \{ai + a^y ui - ai 

= ^/^{crl + a3)U{ai,a3) ( 6 . 11 ) 


Thus we find a string equation for three marked points of the Kontsevich- 
Penner model, 


^ ^ Pii —1P12 ^g T ^ ^ni^n2 —1 >g . (6.12) 


Repeating the same procedure, we have a string equation of s-marked 
points for the Kontsevich-Penner model, 

S S S 

< ^0 n >9= < ^n,-l n >9 (6-13) 

i=l j=l i^j 


W-constraints equation 


We consider next the terms of order a 2 in the expression 2 marked points 
U(ai,a 2 ). The term <72 corresponds to ti ~ tr^. Such fractional indices 
correspond to W constraints [19], which appear in the p-th higher Airy matrix 
model. Such a fractional index appears also in the non-orientable Lie algebra 
0{2N) as we have seen. Therefore, the terms t^_^i are characteristics of open 
intersection numbers. Keeping the order a 2 terms in (6.5), and putting the 
other (72 to zero. 


K((7i, (72) = (72i/aie J dadp 

{a + (d){iui) 
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duidu 2 -ui-ul+k\ogi - 
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(27ri)2 
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(6.14) 


where we use the following contour integrals, 
r dii. 'I 9 1 

R,2n 

J 'Mil M""" Z 

du 
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'2n+l 


27ri n2n+i 
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p^-2 / 27ri n2n+i 


R 2 — —2Rq — — 2r(-) 


2 

lim -r(l 
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2n 
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(6.15) 
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From this expression, we have 


7=i= k, 


< TiT3 
2 2 


3 >3=1= < riTs >3=3= -(3/c + Ak^) (6.16) 


dilaton equation 

Next, we consider the dilaton equation, which involve as, i.e. ri in the 

2 

intersection numbers. The equation LqZ = 0 is a dilaton equation, which is 
derived by considering the terms of order tr^. For the dilaton equation, we 
have 

S S 

< n n > 9 = ( 2 fi' - 2 + s) < n >g ( 6 . 17 ) 

i=l i=l 

From the three point function U{ai, <72, <73), we obtain 

< ToTin >3=i = < Ton >3=1= k, < ToVi >3=0 = < >g= 0 = 1 , ••• (6.18) 

which satishes the dilaton equation of (6.17). 

3 1 

For (72 in U{ai, (72, ..., (7^), scaling Ui —)■ (—) 2 , it is obvious that the loga¬ 
rithmic term for U2 can be neglected, since it gives higher orders. Therefore, 
the dilaton equation does not show the effect of the logarithmic term and the 
equation is same as the dilaton equation without the logarithmic potential 
(there is no k in the equation) as (6.17), and thus we find that (6.17) holds. 

In Appendix, we discuss the string equation, the divisor equation and the 
dilaton equation for p spin curves in the presence of the logarithmic potential. 

Virasoro equations for Kontsevich-Penner model 

The Virasoro equations are expressed through operators L„, which act 
on the partition function Z = , 

LnZ = Lne^ = 0 (6.19) 


with n = —1, 0,1,... ; for n = —1 it gives the string equation and for n = 0 
the dilaton equation. The free energy is a generating function of the inter¬ 
section numbers with a variables = tr—f-r. For the Kontsevich-Penner 

A“+7 

model, since there is a logarithmic potential, we need to consider also half¬ 
integer values for n in We recall previous calculations [19] and summarize 
here a comparison with the calculation based7/((7i,..., (7^). The first Virasoro 

equation for the order tr^, which is the string equation, is given by 

A 2 


dF 1 2 


k 

— t 1 + 

2 2 
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n 1 




dF 

dtn 


( 6 . 20 ) 
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The free energy is a generating function of the intersection numbers as 


fdn 

d„ n n 


( 6 . 21 ) 


where 


We have the relation 


tr, = tr- 


(23A)^ 


( 6 . 22 ) 


K = (23)"+2 4 (6.23) 

In [19], the Virasoro equations up to third order (i.e. up to the dilaton 
equation) are obtained as 

(-^ + - ^ti)g = 0 {string equation) 


.A 

‘ dti 


1 k'^ 


iA’ + - l/hs = 0 


(-3—-- Tr ~ ^toti — — -J^^l)g = 0 {dilaton equation) 


(9ti 


16 


(6.24) 


with 


J« = 
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mx-rn, (m = ...,-2,-1,0,1,2,...) 


dXm 

j{2) ^ y- . t(1) t{1) 

= E w^ + 2 E E 


i+j=m 9XidXj 


tW t(1) t(1) . 


‘ dx-i 




7(3) _ . t( 1 ) t( 1 ) t( 

., dxidxjdxh ^ 

i+j+k=m ^ J 


(6.25) 


—i-\-j-\-k=—4 

where : ■ ■ ■ : means normal ordering, i.e. pulling the differential operator 
to the right. Xn = The solution of the Virasoro equations, which 

includes the string equation, gives [19] 


F — —tl + —{l + 12k‘^)ti + -ktoti 


(6.26) 


+ ^^0^1 + + ^(1 + 12A;2)tot2 

k 1 1 

H— t^ts T —/unfits T —{k T k^)t6 T ■ ■ ■ 

4^2 4 22 6 ^ ^2 


(6.27) 
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From this expression, the intersection numbers , which are dehned as 
(6.21) are obtained by changing tn to tn for small genus, 

1 + 12k'^ 

< >5=0= 1, < n >g=i= — — —, < Tori >g=i= k, 

1 + 12A;^ 

< ToVi >g= 0 = 1 , < >3=1= - — -, < ToTiTi >3=1= k, 

1 + 12/c^ 

<ror2>3=i =—— -, <T^Tl>g^l=k, 

< TITS >3=1= /c^, < rs > 0 - 3 = -^(/c + k^) ( 6 . 28 ) 

These intersection numbers satisfy the string equation, the W-constraints, 
and the dilaton equation. They provide identical values as those calculated 
from U{ai, ...,0's) for the Kontsevich-Penner model, as shown in Appendix. 

7 Gromov-Witten invariants of CP^ model 


The Gromov-Witten invariants of CP^ model has been studied [46, 47, 48]. 
Recently, the Gromov-Witten invariants are evaluated in more higher orders 
[49]. We apply the present method to Gromov-Witten invariants of CP^ 
model, since it has a similar matrix model representation as Kontsevich type 
of the external source [48]. 

The CP^ matrix model is described as 


^ — J ^)+trMA 


We use the Gaussian random matrix model with an external source, and by 
the tuning the external source A, we obtain (7.1) as a generalized Kontsevich 
model. Therefore, as before, we consider the Fourier transform of the density 
correlation functions. Particularly, we consider U{a), which is 


a I 2Tii 




1 f du 


,(2sinh|)(e“-qe-“)+u 


a J 2m 


By a; = e“, we have 


U{a) = - (f -^g(2>3inh^)G-5a; 1) 

a J 2ni 


where we inserted N to make clear the genus expansion. The residue calcu¬ 
lation becomes 

(7,4) 
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In the genus zero, N —)■ oo, dropping the irrelevant factor of q, we have 


UM = E 




a 


2{d-l) 


= E M = E < >5=0 (d + l)a 


2d 


d^o (d + 1)W! 


d=0 


with 


< T2d >9=0— 


(d+l)!2' 


(7.5) 


(7.6) 


For the higher genus g, we expand (2iVsinh^)^“ \ and pick up the genus 
g terms from order terms. 


-(2iVsinh— 2 d^ i‘^d l)(10d 7) 

2iV"" tr^criATd ' ' 


24iV2 


5760iV4 


With the shift of power of a, 2d + 2n 2d, we obtain Gromov-Witten 
invariants of genus g as 


with 


t'(<T) = EE < T2d >g {d + 1 — g)a 
d 9 


2d 


(7.8) 


< r2d{t^) >3=0 = 

< r2d{u}) >3=1 = 

< T2d(w) >3=2 = 

< T2d{uj) >3=3 = 


((d+l)!)2 
2d - 1 


24(d!)2 

d2(2d-3)(10d- 17) 

27.32.5(d!)2 

d2(d - l)2(2d - 5)(140d2 - 784d + 1101) 
210 ■ 34 • 5 ■ 7(d!)2 


(7.9) 


These numbers agree with the result of the recent evaluation by Norbury 
and Scott by a different method up to genus three[49]. It is straight forward 
to evaluate Gromov-Witten one point invariants in any order of genus from 
U{a). We have for = 4 as 


< r2d{uj) >3=4= 


d^{d - lf{d - 2)2 (2d - 7)(10d - 39)(140d2 - 1092d + 2143) 


(d!)^ 


215 • 35 • 52 ■ 7 


(7.10) 
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8 Discussions 


In this article, we have considered the generalization of the Airy matrix 
model to a p-th singnlarity. This provides the intersection nnmbers of the 
moduli space of p-spin curves for orientable and non-orientable Riemann 
surfaces, with Lie algebras of 17(A^), 0{2N), 0(2A^+1) and Sp{N). The Euler 
characteristics are easily evaluated by taking the p —)■ — 1 limit. Our results 
are consistent with the two categories, orientable and non-orientable surfaces, 
since we have obtained two type of topological invariants (two different Euler 
characteristics) for Lie algebras. The expressions agree with the virtual Euler 
characteristics obtained earlier [42] for non-orientable surfaces. We have 
obtained explicit expressions to all order in the genus for one marked point 
in the p = 3 and p = 4 cases given in terms of Bessel functions. 

For the open intersection numbers, which are dehned by the insertion 
of a disk on a closed Riemann surface as a boundary, we have used the 
Kontsevich-Penner model. We have derived the Virasoro equations, string 
equation and dilaton equations, for this Kontsevich-Penner model from ex¬ 
plicit integral representations. The open intersection numbers are extended 
to p-spin curves, from a higher Airy matrix model with a logarithmic poten¬ 
tials. 

In our previous article [6], the Airy matrix model with a logarithmic 
potential was derived from the average of two characteristic polynomials in a 
two matrix model with an external source The eigenvalues of the first matrix 
Ml is on an edge of the distribution, and for the other matrix M 2 in the bulk,. 
After integration over the matrix M 2 , a model with a logarithmic potential is 
obtained for Mi. The coefficient of the logarithmic potential k corresponds 
to the power of (detMi)^. This logarithmic potential provides the boundary 
for the open intersection theory. 

The integral representation of the s point correlation function for a Gaus¬ 
sian matrix model with an external source provides a powerful tool for the 
evaluation of open/close intersection numbers and Gromov-Witten invari¬ 
ants. It would be interesting to extend the present analysis to more compli¬ 
cated cases , such as the Gromov-Witten theory of CP^~^. 
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Appendix: Virasoro equations of open intersection 
numbers for p spin curves 

The open/close intersection numbers of p-spin curves are evaluated from 
the expressions for the s-point correlation functions U{ai, ...,as)- From 
(5.26), the one marked point intersection numbers < >g are obtained 
as 


< n,o > 1 = 


P-1 + I2k‘^ 
24 


< Tnj >1= —{pk + 2k^) 
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t 24'"" ' " 'r(i - 


^ '^n,j ^2 — 
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< Tnj >5= {-(18 — 25p + 30p^ — 5p^) H- ip — 1)}- 

2 ^5760^ ^ ^ 144 

< TnJ >3= 4i onnL.n (P “ 1)(P “ 5)(1 + 2p){8p‘^ - 13p - 13) 


X- 


p2'2903040 

— ^ —(—lOp^ + 85p^ + 90p + 19)/c^ H-—(5p + 3)k* H- ^ — k^} 

57600^ ^ ^ y J 3gQQ J 

r(i-!) 


r(i - 


i+i' 


where n and j are constrained by the condition, 

(p+ l)(2p- 1) =pn + i + l 


(A.l) 

(A.2) 


i.e. he s = 1 case of the general condition, 

S S 

(P+l)(2fi'- 2 + s) = ^ jfc + s (A.3) 

k=l k=l 

For < rs >„=3 of p = 2 is obtained from lim „^2 < 'r 2 i >„=3 since the right 
hand side of (A.2) is the same, and it becomes ^{k + k^)- 

For the 2 marked points, the intersection numbers of the p spin curves 
are derived from U{ai,a 2 ). 


f/(ai,a2) = 


/ duidu2 + 1 ( 7 , 


J {27ii) 


n 


_i Ui rfC^i 


i=l 
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(til -U2 + \{(Ji + (T2)){u2 - Ml + \{(Ji + 0 - 2 )) 


(A.4) 


string equation 


Using Xi = aiuf, and a —)■ {<Jia 2 )pa, j3 —)■ (ai(T 2 )^/d, taking the same 
process as (6.5), we obtain 


U{ai,a2) = 


. P(P-I) 


e 24 
e 


JO 

2 +— 1 —— 

E.Ui ")+■■■ 

h 


dxidx2dadj3{xiX2)p e 


— 1 ^ — X^ — 


X 1 —X 2 


ict((CT2tci)f'—(f 7 iX 2 )P+|(o-if 72 )P (fTi+fT 2 ))g^/ 9 ((fTia; 2 )f'—(fT 23 :i)P+1(0-10-2)P (0-1+0-2)) 

^ 1+i 

2 x? + i(T- 

' 2 ^ \k 

*=' + - ia, 


1+1 ■ 


(A.5) 


Taking the term of order a 2 and neglecting higher order terms in < 72 , we have 

U{ai, a2) = r(l --)■</' 

~ v J 27 ri Ml — l(Ti 


TT 


= -^r(l--)-aiU(ai) 

TT p 

This equation is a string equation, 


(A.6) 


(A.7) 


A string equation for three marked point for p spin curves is an extension of 
(6.11). It is easily obtained from 


U{ai,a 2 ,a 3 ) = a2'’r(l - -)(ai + a3)U(ai, 0-3) 

p 


(A.8) 


which is the string equation, 

AtjO+nijil-n2,j2 ^g ^ Ail—i,jiTi2,i2 ^g ^ Aii,iiAi2— 1,12 ^g (A-9) 


W constraint equation 

We consider next the W constraint equation. Since there are spin j = 

0,1, ■■■,p — 1 indices for the intersection numbers of p spin curves, we have an 

2 

equation which involves tqj. Taking next al, and neglecting higher terms in 
(72 in (A.5), we obtain the intersection number with Toy. 
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The two point correlation function U{ai,a 2 ) is expressed as [7], by using 
the following representation. 

1 1 
U1-U2 + \{ai + (T2) U1-U2- |(cTi + (T2) 

=-i- r dae-“(“i-“ 2 )sinh(^(ai + ^2)) (A.IO) 

(Tl+ (J2 Jo 2 

By a -)■ ((Ti(T2)pa, Ui = (xi/a*)?, sinh(f (ai + (J2)) ~ f (ai0-2)5(ai + 0-2), we 
1 

obtain af term as 


77 ( 0 - 1 , 0 - 2 ) = 


1 ( 0 - 10 - 2 ) p 




dxidx2- 


( 0 - 2 X 1 ) 


(X 1 X 2 ) 




= o-f 0-277(0-2) 


(A. 11 ) 


which is a string equation. For af term, we expand e = 1 + (o-iX 2 )j’Q;. 

77 (ai,a 2 ) =-afr(l—)a 2 M dx 2 x| e 24 ^"2 ^2 + (_- 

p p Jo ^p _l 


X 2 - 2^2 


(A.12) 


From above integral, we obtain o-f 0-2 ^ term as 


jT(-r ft rrn 2 ip(p-i)2 (p-i)(p-2)(p-3)^ 

£/(<Ti.<T2) ~ « 0 ,ii 4 .i|r(l - -)1 {(1 - -5P-} 

(A.13) 

Separating a factor (1 — ^), we obtain in the case k = 0, 


< 'TopTi,! >g — < '^3,2 >g A Ti^q >g 

for general p. Similarly we obtain for g = 3, when k = 0, 

1 

< Xo, 1 X 6,3 >g = < X5,4 >g=3 +” < Xi,o >g=l< T3,2 >g=2 


(A.14) 


(A.15) 


dilaton equation 

The dilaton equation for p-spin curves is 


A B,0 n Fifcjfc {‘J‘9 2 + s) < J([ Tn^.j^ >g 


k=l 


k=l 


(A.16) 
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We consider s = 1, two point correlation function U{( 71 ,( 72 ). By the shift 
a —)■ (ai(T 2 )p, Xi = 7iu\, we have 


U { 71 , 72 ) = 


7i + 72 


n 


POO i i 2-1- 

/ dasinh(-(aia2)2(ai+ ^ 

JO 2 

(A.17) 


3 

.Xi^\7j 

3 J 

Xi - ^7l 


For simplicity, we evaluate the p = 2 case. The term of order 7i comes from 

1 O' . s i 1 3 1 3 7 

sinh(-(aia2)(ai + o-| (A.18) 


7i + (72 


By the integrations of a and Xi, we obtain the order of 7 1 

U{71,72) ~ 2afr(l - ^)[72 " J da;2x|e""^2+'" + ^7^ J 6x2X2 

(A.19) 

Noting the integral by parts for the hrst term, we have 


17 { 7 i , 72 ) = 2afr(l - ^)(1 + ^ 72 ) 1 /{ 72 ) 


Since we have 


C7(ff2) = E 


2 (12)’9! 

we resum the two terms of (A.20) as 


-1)^^2 


3ff-i 


TT(rr rr 1 = 9n-^ Fl'l'l ( 1)^ 3(g-|) ( 1)^^^ 1 3 (g+i)^ 

( 1’ 2 ) ^ (12)J^! 6 ^ 

which provides the dilaton equation for p = 2, 

^ B oFfjj' {2.^ 1) ^ Aij ^g 

We can check , for instance, g = I case for p = 2 as 

< rjo >g=i=< ri,o >g=i= ^(1 + I2k^) 


(A.20) 

(A.21) 

(A.22) 

(A.23) 

(A.24) 


The above equation may be easily extended to p > 2 by the same process. 
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